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Abstract — We use information theoretic achievable rate formu- 
las for the multi-relay channel to study the problem of optimal 
placement of relay nodes along the straight line joining a source 
node and a destination node. The achievable rate formulas that 
we utilize are for full-duplex radios at the relays and decode- 
and-forward relaying. For the single relay case, and individual 
power constraints at the source node and the relay node, we 
provide explicit formulas for the optimal relay location and the 
optimal power allocation to the source-relay channel, for the 
exponential and the power-law path-loss channel models. For the 
multiple relay case, we consider exponential path-loss and a total 
power constraint over the source and the relays, and derive an 
optimization problem, the solution of which provides the optimal 
relay locations. Numerical results suggest that at low attenuation 
the relays are mostly clustered close to the source in order to be 
able to cooperate among themselves, whereas at high attenuation 
they are uniformly placed and work as repeaters. We also prove 
that a constant rate independent of the attenuation in the network 
can be achieved by placing a large enough number of relay nodes 
uniformly between the source and the destination, under the 
exponential path-loss model with total power constraint. 

I. Introduction 

Wireless interconnection of mobile user devices (such as 
smart phones or mobile computers) or wireless sensors to 
the wireline communication infrastructure is an important 
requirement. These are battery operated, resource constrained 
devices. Hence, due to the physical placement of these devices, 
or due to the channel conditions, a direct one-hop link to 
the infrastructure "base-station" might not be feasible. In 
such situations, other nodes could serve as relays in order 
to realize a multi-hop path between the source device and 
the infrastructure. In the cellular context, these relays might 
themselves be other users' devices. In the wireless sensor 
network context, the relays could be other wireless sensors 
or battery operated radio routers deployed specifically as 
relays. In either case, the relays are also resource constrained 
and a cost might be involved in engaging or placing them. 
Hence, there arises the problem of optimal relay placement. 
Such an optimal relay placement problem involves the joint 
optimization of node placement and of the operation of the 
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Fig. 1 . A source and a destination connected by a multi-hop path comprising 
N relay nodes along a line. 

resulting network, where by "operation" we mean activities 
such as transmission scheduling, power allocation, and channel 
coding. 

In this paper we consider the problem of maximizing the 
data rate between a source node (e.g., a wireless sensor) and a 
destination (e.g., a sink node in a wireless sensor network) by 
means of optimally placing relay nodes on the line segment 
joining the source and the destination; see Figure [T] In order 
to understand the fundamental trade-offs involved in such a 
problem, we consider an information theoretic model. For a 
placement of the relay nodes along the line and allocation of 
transmission powers to these relays, we model the "quality" 
of communication between the source and the destination by 
the information theoretic achievable rate of the relay channeQ 
The relays are equipped with full-duplex radio^ and carry out 
decode-and-forward relaying. We consider scalar, memoryless, 
additive white Gaussian noise channels. , and carry out decode- 
and-forward relaying. We consider scalar, memoryless, time- 
invariant, additive white Gaussian noise (AWGN) channels. A 
path-loss model is important for our study, and we consider 
both power-law and exponential path-loss models; the channels 
are assumed to be time invariant. 

A. Related Work 

A formulation of the problem of relay placement requires 
a model of the wireless network at the physical (PHY) 
and medium access control (MAC) layers. Most researchers 
have adopted the link scheduling and interference model, i.e., 
a scheduling algorithm determines radio resource allocation 
(channel and power) and interference is treated as noise 
(see |3|). But node placement for throughput maximization 
with this model seems to be intractable because the optimal 

'The term "relay channel" will, in this paper, include the term "multi-relay 
channel." 

^See (T|, f2l for recent efforts to realize practical full-duplex radios. 



throughput is obtained by first solving for the optimum sched- 
ule assuming fixed node locations, followed by an optimization 
over those locations. Hence, with such a model, there appears 
to be little work on the problem of jointly optimizing the 
relay node placement and the transmission schedule. In (4\, 
the authors considered placing a set of nodes in an existing 
network such that certain network utility (e.g., total transmit 
power) is optimized subject to a set of linear constraints 
on link rates. They posed the problem as one of geometric 
programming assuming exponential path-loss, and showed that 
it can be solved in a distributed fashion. To the best of our 
knowledge, there appears to be no other work which considers 
joint optimization of link scheduUng and node placement using 
link scheduling model. 

On the other hand, an information theoretic model for a 
wireless network often provides a closed-form expression for 
the channel capacity, or at least an achievable rate region. 
These results are asymptotic, and make idealized assumptions 
such as full-duplex radios, perfect interference cancellation, 
etc., but provide algebraic expressions that can be used to 
formulate tractable optimization problems. The results from 
these formulations can provide useful insights. In the context 
of optimal relay placement, some researchers have already 
exploited this approach. For example, Thakur et al. in fSl 
report on the problem of placing a single relay node to 
maximize the capacity of a broadcast relay channel in a 
wideband regime. The linear deterministic channel model (|6|) 
is used in [7] to study the problem of placing two or more 
relay nodes along a line so as to maximize the end-to-end 
data rate. Our present paper is in a similar spirit; however, we 
use the achievable rate formulas for the A^-relay channel (with 
decode and forward relays) to study the problem of placing 
relays on a line under individual node power constraints as 
well as with sum power constraints over the source and the 



single relay placement problem in this context (Theo- 
rem |5]). A numerical study shows that, the relay nodes 
are clustered near the source at low attenuation and are 
placed uniformly between the source and the destination 
at high attenuation. We have also studied the asymptotic 
behaviour of the achievable rate i?^r when N relay nodes 
are placed uniformly on a line of fixed length, and show 
that for a total power constraint Pt among the the source 
and the relays, liminfjv-j-oo ^^jv > where C(-) 

is the AWGN capacity formula and cr^ is the power of 
the additive white Gaussian noise at each node. 

The rest of the paper is organized as follows. In Section 
n] we describe our system model and notation. In Section 
Tn) node placement with per-node power constraint has been 
discussed. Node placement for total power constraint has been 
discussed in Section HV] Conclusions have been drawn in 
Section |V] The proofs of all the theorems are given in the 
appendices. 



II. System Model and Notation 



The multi-relay channel was studied in (S) and JU and is 
an extension of the single relay model presented in ifTOl . We 
consider a network deployed on a line with a source node, a 
destination node at the end of the line, and N full-duplex relay 
nodes as shown in Figure [T] The relay nodes are numbered 
as 1,2, •• • ,N. The source and destination are indexed by 
and iV + 1, respectively. The distance of the k-th node from 
the source is denoted hy yk := J'l + ^2 + • • • + r^. Thus, 
Un+i = L. As in ||8l and ||9], we consider the scalar, time- 
invariant, memoryless, additive white Gaussian noise setting. 
A symbol transmitted by node i is received at node j after 
multiplication by the (positive, real valued) channel gain hij. 
The Gaussian additive noise at any receiver is independent and 
identically distributed from symbol to symbol and has variance 
In Section III we consider the problem of placing a single a^. The power gain from Node i to Node j is denoted by 

9i,j ~ We model the power gain via two alternative path- 



relays. 



B. Our Contribution 



relay with individual power constraints at the source and 
the relay. In this context, we provide explicit formulas for 
the optimal relay location and the optimal source power 
split (between providing new information to the relay and 
cooperating with the relay to assist the destination), for 
the exponential path-loss model (Theorem |2]i and for the 
power law path-loss model (Theorem [3]). We find that at 
low attenuation it is better to place the relay near the 
source, whereas at very high attenuation the relay should 
be placed at half-distance between the source and the 
destination. 

In Section IV we focus on the N relay placement 
problem with exponential path-loss model and a sum 
power constraint among the source and the relays. For 
given relay locations, the optimal power split among the 
nodes and the achievable rate are given in Theorem |4] 
in terms of the channel gains. We explicitly solve the 



loss models: exponential path-loss and power-law path-loss. 
The power gain at a distance r is e^'"' for the exponential path- 
loss model and for the power-law path-loss model, where 
p > 0, ?7 > 1. These path- loss models have their roots in the 
physics of wireless channels (ifTTl). For the exponential path- 
loss model, we will denote A pL. Under the exponential 
path-loss model, the channel gains and power gains in the line 
network become multiplicative, e.g., ft.i,i+2 — /ii.i+i/ii4.i,i+2 
and 5i^j+2 = g<,j+i5i+i,»+2 for i e {0, 1, • • • ,N - 1}. In this 
case, we define i = 1 and hi^i — 1. The power-law path-loss 
expression fails to characterize near-field transmission, since it 
goes to oo as r — > 0. One alternative is the "modified power- 
law path-loss" model where the path-loss is minjr^'', 6^^'} 
with > a reference distance. In this paper we consider both 
power-law and modified power-law path loss models, apart 
from the exponential path-loss model. 



A. The Multi-Relay Channel 

For the multi-relay channel, we denote the symbol trans- 
mitted by the i-th node at time t (t is discrete) by Xi{t) 
for i = 0, 1, • • • , TV. Zk{t) ^ 7V(0, cr^) is the additive white 
Gaussian noise at node k and time t, and is assumed to be 
independent and identically distributed across k and t. Thus, 
at symbol time t, node fc,l <A:<A^+1 receives: 

Yk{t)= J2 hj^kXjit) + Zk{t) 

ie{o,i,--- ,N}.j^k 

hi m, the authors showed that an inner bound to the capacity 
of this network is given by (defining C{x) ^ log2(l + x)): 

^ ^ ^ j = l i=0 ' 

where Pi ^ denotes the power at which node i transmits to 
node j. Its significance will be clear later in this section. 

In order to provide insight into the expression in Equation 
([TJ and the relay placement results in this paper, we provide 
a descriptive overview of the coding and decoding scheme 
described in ISj. Transmissions take place via block codes of 
T symbols each. The transmission blocks at the source and the 
iV relays are synchronized. The coding and decoding scheme 
is such that a message generated at the source at the beginning 
of block 6, 6 > 1, is decoded by the destination at the end of 
block b + N, i.e., 1 block durations after the message was 
generated (with probability tending to 1, as T — > oo). Thus, 
at the end of B blocks, _B > iV + 1, the destination is able to 
decode B ~ N messages. It follows, by taking B — > oo, that, 
if the code rate is R bits per symbol, then an information rate 
of R bits per symbol can be achieved from the source to the 
destination. 

As mentioned earlier, we index the source by 0, the relays 
by fc, 1 < k < N, and the destination by + 1. There 
are {N + 1)^ independent Gaussian random codebooks, each 
containing 2^^^ codes, each code being of length T; these 
codebooks are available to all nodes. At the beginning of 
block 6, the source generates a new message Wb, and, in this 
description, we assume that each node fc, 1 < k < N + 1, has 
a reliable estimate of all the messages w^^j,] > fc. In block 
h, the source uses a new codebook to encode w^. In addition, 
relay k,l < k < N, and all the transmitters previous to it 
(indexed < j < k — 1), use another codebook to encode 
Wb-k (or their estimate of it). Thus, if the relays 1, 2, • • • ,k 
have a perfect estimate of Wb-k at the beginning of block b, 
they will transmit the same codeword for Wb-k- Therefore, in 
block 6, the source and relays 1, 2, • • • , /c coherently transmit 
the codeword for Wb-k- In this manner, in block b, transmitter 
k,0 < k < N, generates + 1 — fc codewords, corresponding 
to Wb-kjWb-k-i, ■ ■ ■ , Wb^N, which are transmitted with pow- 
ers Pk,k+i,Pk,k+2, ■■■ , Pk,N+i- In block b, node fc, 1 < fc < 

+ receives a superposition of transmissions from all other 
nodes. Assuming that node fc knows all the powers, and all the 



channel gains, and recalling that it has a reliable estimate of 
all the messages Wb-j,j > fc, it can subtract the interference 
from transmitters fc + 1, fc + 2, • • • , A^. At the end of block b, 
after subtracting the signals it knows, node fc is left with the 
fc received signals from nodes 0, 1, • • • , (fc — 1) (received in 
blocks 6, 6 — 1, • • • ,b — k + 1), which all carry an encoding 
of the message Wb-k+i- These fc signals are then jointly 
used to decode Wb-k+i, using joint typicality decoding. The 
codebooks are cycled through in a manner so that in any block 
all nodes encoding a message (or their estimate of it) use the 
same codebook, but different (thus, independent) codebooks 
are used for different messages. Under this encoding and 
decoding scheme, relatively simple arguments lead to the 
conclusion that any rate strictly less than R displayed in 
Equation ([T} is achievable. 

From the above description we see that a node receives 
information about a message in two ways (i) by the message 
being directed to it cooperatively by all the previous nodes, 
and (ii) by overhearing previous transmissions of the message 
to the previous nodes. The higher the channel attenuation, 
the less will be the contribution of farther nodes, "overheard" 
transmissions become less relevant, and coherent transmission 
reduces to a simple transmission from the previous relay. The 
system is then closer to simple store-and-forward relaying. 

The authors of ^ have shown that any rate strictly less 
than R is achievable through this coding and decoding scheme 
which involves coherent multi-stage relaying and interference 
subtraction. This achievable rate formula can also be obtained 
from the capacity formula of a physically degraded multi-relay 
channel (see f9l for the channel model), since the capacity 
of the degraded relay channel is a lower bound to the actual 
channel capacity. In this paper, we will seek to optimize R over 
power allocations to the nodes and the node locations. Also, 
log(-) in this paper will mean the natural logarithm unless the 
base is specified. We denote the value of R optimized over 
power allocation and relay locations by R*. 

For the single relay channel, A^ = 1. Thus, by ([T]), an 
achievable rate is given by (see also IJOJ ): 

D ■ / r< f 9o.iPo,i\ 

The following result justifies our aim to seek optimal relay 
placement on the line joining the source and the destination 
(rather than anywhere else on the plane). 

Theorem 1: For given source and destination locations, the 
relays should always be placed on the line segment joining the 
source and the destination in order to maximize the end-to-end 
data rate between the source and the destination. 

Proof: See Appendix [A] ■ 
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III. Single Relay Node Placement: Node Power 
Constraints 

In this section, we aim at placing a single relay node 
between the source and the destination in order to maximize 
the achievable rate R given by Equation Q. Let the distance 
between the source and the relay be r, i.e., ri ~ r. Let 
a := We assume that the source and the relay use the 
same transmit power P. Hence, Pq.i = o^Pj Pq,2 = (1 — ct)P 
and Pi 2 — P- Thus, for a given placement of the relay node, 
we obtain: 



R = min 



C 



V <t2 



C(^^ (^30,2 + 31,2 + 2^(1 - 0)90,291,2^^ I (5) 

Maximizing over a and r, and exploiting the monotonicity 
of C( ) yields the following problem: 



max max 

relO.L] aG[0 



c^^ min |ago,i, 90,2 + 9i,2 + 2^ (1 - 0)90, 291, 2| 



(6) 



Note that 50,1 ™d gi 2 in the above equation depend on r, 
but 30,2 does not depend on r. 
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A. Exponential Path Loss Model 

1) Optimum Relay Location: Here (70,1 = 
g-piL-r) ^jjjj _ e^P^. Let A := pL, and the optimum 
relay location be r* . Let x* := ^ be the normalized optimal 
relay location. Let a* be the optimum value of a when relay 
is optimally placed. 

Theorem 2: For = 1, under the exponential path-loss 
model, if the source and relay have same power constraint 
P, then there is a unique optimum relay location r* and the 
following hold: 

(i) X* — max{a::+, 0} where x+ — log ^2e~''' + e^^^ 

(ii) For < A < log 2, x* ^0,a* = 1 and R* = C(^). 

(iii) For log2 < A < log 4, x* = 0, a* = 4e"^(l - e"-^) and 



R* 



C(-^4e-^(l-e-^)). 



(iv) For A > log 4, x* ^ x^. Moreover, a* and R* are given 
by Equations ([3]) and Q. 

Proof: See Appendix [B] ■ 

2) Numerical Work: In Figure |2] we plot x* and a* 
provided by Theorem |2] versus A. Recalling the discussion 
of the coding/decoding scheme in Section |llj we note that the 
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Fig. 2. Single Relay, exponential path-loss, node power constraint: x* and 
a* versus A. 

relay provides two benefits to the source: (i) being nearer to 
the source than the destination, the power required to transmit 
a message to the relay is less than that to the destination, and 
(ii) having received the message, the relay coherently assists 
the source in resolving ambiguity at the destination. These two 
benefits correspond to the two terms in Equation (|2]). Hence, 
Theorem [2] and Figure [2] provide the following insights: 

• At very low attenuation (A < log 2), x* — (since x^- < 
0). Since gi 2 > 5 and 50.2 > 5, we will always have 
50,2+31,2 + 2^/(1 - a)go,2gi,2 > ago,i for all a € [0, 1] 
and for all r € [0,L]. The minimum of the two terms in 
Equation (j6|l is maximized by a = 1 and r — Q. 

• For log 2 < A < log 4, we again have x* < 0. a* 
decreases from 1 to | as A increases. In this case since 
attenuation is higher, the source needs to direct some 
power to the destination for transmitting coherently with 
the relay to the destination to balance the source to 
relay data rate. A = log 4 is the critical value of A at 
which source to relay channel ceases to be the bottleneck. 
Hence, for A € [log 2, log 4], we still place the relay at 
the source but the source reserves some of its power 
to transmit to the destination coherently with the relay's 
transmission to the destination. As attenuation increases, 
the source has to direct more power to the destination, 
and so a* decreases with A. 

• For A > log 4, X* > 0. Since attenuation is high, it is 
no longer optimal to keep the relay close to the source. 
Thus, X* increases with A. As A — > 00, a;* — > ^ and 
a* — 7> 1. We observe that the ratio of powers received 
by the destination from the source and the relay is less 
than e^2, since x* < ^. This ratio tends to zero as A — > 
00. Thus, at high attenuation, the source transmits at full 
power to the relay and the relay acts just as a repeater 

• We observe that x* < ^. The two data rates in Equation 
(jsjl can be equal only if r < ^ . Otherwise we will readily 
have (70,1 ^ 9o,2 + .91.2, which means that the two rates 
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Fig. 3. Power law path-loss, single relay node, individual power constraint: 
X* and a* versus r]. 



Fig. 4. Modified power-law path-loss (with -j- = 0.1), single relay node, 
individual power constraint: x* and a* versus rj. 



will not be equal. 



IV. Multiple Relay Placement : Sum Power 
Constraint 



B. Power Law Path Loss Model 



In this section, we derive the optimal placement of relay 
nodes to maximize R (see Equation ([T]l), subject to a total 
power constraint on the source and relay nodes given by 



1) Optimum Relay Location: For the power-law path-loss X^i^o ~ ^t- We consider only the exponential path 



go.2 = i and gi,2 = [L — r) Let x* loss model. We will first maximize R in Equation ([T]) over 

Pj J , < i < j < {N + 1) for any given placement of nodes 



model, (7o,i = 

be the normalized optimum relay location which maximizes 
R. Then the following theorem states how to compute x*. 

Theorem 3: (i) For the single relay channel and power- 
law path-loss model (with rj > 1), there is a unique place- 
ment point for the relay on the line joining the source and 
the destination, maximizing R. The normalized distance 
of the point from source node is x* < ^, where x* 
is precisely the unique real root p of the Equation 

_ ^^2(1 _ (1 _ 1)-,) ^ (1 _ _ (1 _ 

in the interval (0, |]. 

(ii) For the "modified power-law path-loss" model with 
2b < L, the normalized optimum relay location is 

X* = max{p, j}. 

Proof: See Appendix [B] ■ 

2) Numerical Work: The variation of x* and a* as a 
function of rj for the power-law path loss model are shown 
in Figure [5] As 77 increases, both x* and a* increase. For 
large rj, they are close to 0.5 and 1 respectively, which means 
that the relay works just as a repeater At low attenuation the 
behaviour is different from exponential path-loss because in 
that case the two rates in Equation (j5]l can be equalized since 
channels gains are unbounded. For small 77, relay is placed 
close to the source, go,i is high and hence a small a suffices 
to equahze the two rates. Thus, we are in a situation similar 
to the A > log 4 case for exponential path-loss model. Hence, 
a* and x* increase with A. 

On the other hand, the variation of x* and a* with 77 
assuming the modified power-law path-loss model are shown 
in Figure|4] with = 0.1. Here, for small values of 77, p < 0.1 
and hence x* ~ 0.1. Beyond the point where 0.1 is the solution 
of (x-^+i - 1)2(1 - (i - l)--)) = (1 - x)-^ - (i - 1)-", 
the behaviour is similar to the power-law path-loss model. 
However, for those values of 77 which result in p < j^, the 
value of a* decreases with 77. This happens because in this 
region, if 77 increases, 50,1 remains fixed but 50,2 decreases. 



(i.e., given yi,y2,''- jVn)- This will provide an expression 
of achievable rate in terms of channel gains, which has to 
be maximized over 7/1, 772, • • • , Z/at- Let 7fc := Y.^~q Pi,k for 
k E {1, 2, • • • , iV + 1}. Hence, the sum power constraint 
becomes J2k=i 7fc = Pt- 

Theorem 4: (i) For fixed location of relay nodes, the 
optimal power allocation that maximizes the achievable 
rate for the sum power constraint is given by: 




yo<i<j<{N + i) 

if j <i 



(7) 



where 



71 = 



Pt 



1 + 50,1 Lfc=2 



(go,fc-i— go.fc) 
fc-i 1 

90, i 



(8) 



50,1- 



Efc— 1 

(go,3-i-3o,j) 



^1 = g„i 



1 I „ v^w+1 (go,fc-i-go,fc) 
-L + 50,1 2^k=2 „ 



-Pt Vj>2 



1=0 



(ii) The achievable rate optimized over the power allocation 
for a given placement of nodes is given by: 



[yi,y2, 



,yN) = c 



1 

30,1 



-E 



N + l (9o,k-l-9o,k) 
k — 2 1 
90,fc90,fc-l 1^1 = J7r, 



(9) 



Proof: See Appendix [C] I 
Remarks and Discussion: 

• We find that in order to maximize R°p^{yi, 7/2, •• • ,Vn) 
we need to place the relay nodes such that ^ 



(a) N = 2 (b) Af = 3 (c) Af = 5 

Fig. 5. A*^ relays, exponential path-loss model: depiction of the optimal relay positions for TV = 2, 3, 5 for vaiious values of A. 
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Fig. 6. Single relay placement, total power constraint, exponential path-loss: 



and optimum versus 
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is minimized. 



2 ffO,fcSO,fc-l Ef=0^ g,, , 

• From Equation (j?]), it follows that any node k < j will 
transmit at a higher power to node j, compared to any 
node preceding node k. 

• Note that we have derived Theorem |4] using the fact that 
go,fe is nonincreasing in k. If there exists some fc > 1 such 
that 5o.fc = So.k+i, i-e, if k-th and (fc + l)-st nodes are 
placed at the same position, then jk+i = 0, i.e., the nodes 
I < fc do not direct any power specifically to relay k + 1. 
However, relay fc + 1 can decode the symbols received 
at relay fc, and those transmitted by relay fc. Then relay 
(fc + 1) can transmit coherently with the nodes ^ < fc to 
improve effective received power in the nodes j > k + 1. 

A. Optimal Placement of a Single Relay Node 

Theorem 5: For the single relay node placement problem 
with sum power constraint and exponential path-loss model, 
the normalized optimum relay location power allocation 
and optimized achievable rate are given as follows : 

(i) For A < logs, ^ = 0, Po,i = ^^0,2 = Pi,2 = 



e^ + l 2 



\^ and R* 



C 



2Pt 



(ii) For X > log 3, ^ 



Po,2 — 
C 



1 Pi 
1 Ft 



ilog(^ 

Pl,2 = 



+ 1 -1 

_ ^/e^ + l-l p. 



-f0,l - -2-' 

and R* = 



Proof: See Appendix [C] ■ 
Remarks and Discussion: 

• It is easy to check that R* obtained in Theorem |5] is 
strictly greater than the AWGN capacity C (^e~-^) for 
all A > 0. This happens because the source and relay 



transmit coherently to the destination. R* becomes equal 
to the AWGN capacity only at A = 0. At A = 0, we do 
not use the relay since the destination can decode any 
message that the relay is able to decode. 
• The variation of ^ and with A has been shown 

n 

in Figure 6 We observe that limA^oo ^ — \, 
limA-).oo -Po,2 = and limA^o -Po.i = Pt- For large 
values of A, source and relay cooperation provides neg- 
ligible benefit since source to destination attenuation is 
very high. So it is optimal to place the relay at a distance 
^. The relay works as a repeater which forwards data 
received from the source to the destination. 

B. Optimal Relay Placement for a Multi-Relay Channel 
As we discussed earlier, we need to place N relay nodes 



such that 



9o,k - 
Vn < 



go.i 
Vn+i 



N+l 



, _„ ^k-i — 1^ is minimized. Here 



We have the constraint < t/i < ^2 < • • ■ < 
= L. Now, writing = e''^*^, and defining 
zq := 1, we arrive at the following optimization problem: 

N+l 



mm <zi + 



s.t 1 < Zl< 



Zk — Zk-1 



fe=2 



Efc-1 

< ZN < ZN+1 



(10) 



One special feature of the objective function is that it is convex 
in each of the variables zi,Z2, - ■ ■ , z^. The objective function 
is sum of linear fractionals, and constraints are linear. 

Remark: From optimization problem ( [T0| we observe that 
optimum zi, Z2, • • • , zat depend only on A :— pL. Since Zk = 
e'^T, we find that normalized optimal distance of relays from 
the source depend only on A. 

Theorem 6: For fixed p, L and a'^, the optimized achievable 
rate R* for a sum power constraint strictly increases with the 
number of relay nodes. 

Proof: See Appendix [C] ■ 

C. Numerical Work 

We discretize the interval [0, L] and run a search program 
to find normalized optimal relay locations for different values 
of A and N. The result is summarized in Figure |5] We 
observe that at low attenuation (small A), relay nodes are 
clustered near the source node and are often at the source 
node, whereas at high attenuation (large A) they are almost 
uniformly placed along the line. For large A, the effect of 
long distance between any two adjacent nodes dominates the 
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N (number of relays) 



Fig. 7. R* vs N for total power constraint witli = 10. 




N (number of relay nodes) 

Fig. 8. R* vs N for total power constraint with = 100. 



gain obtained by coherent relaying. Hence, it is beneficial to 
minimize the maximum distance between any two adjacent 
nodes and thus multihopping is a better strategy in this case. 
On the other hand, if attenuation is low, the gain obtained 
by coherent transmission is dominant. In order to allow this, 
relays should be able to receive sufficient information from 
their previous nodes. Thus, they tend to be clustered near the 
source 

We plotted R* against the number of relays in Figures [7] and 
|8] The capacity appears to be strictly increasing with number 
of relays, as claimed in Theorem [6] 

D. Uniformly Placed Relays, Large N 

If the relays are uniformly placed, the behaviour of 
R°p^{yiTV2,' ■ ■ ,Vn) (called Rjsj in the next theorem) for 
large number of relays is captured by the following Theorem. 

Theorem 7: For the exponential path-loss model with total 
power constraint, if iV relay nodes are placed uniformly 
between the source and destination, resulting vs\Rj^ achievable 
rate, then liminfTv-Rw > C ^ ^ and limsup^i?Ar < 

cm- 

Proof: See Appendix |C] ■ 

Remark: From Theorem |7] it is clear that we can ensure 
at least C(^^) data rate (roughly) by placing a large enough 
number of relay nodes. 

'At low attenuation, one or more than one relay nodes are placed very 
close to the source. We believe that some of them will indeed be placed at 
the source, but it is not showing up in Figure [s] because we discretized the 
interval [0, L] to run a search program to find optimum relay locations. 



V. Conclusion 

We have studied the problem of placing relay nodes in a line 
network in order to maximize the end-to-end achievable data 
rate. We obtained explicit formulas for the optimal location of 
a single relay node with individual power constraints on the 
source and the relay under power-law and exponential path- 
loss models. For the multi-relay channel with exponential path- 
loss and total power constraint, we derived an expression for 
the achievable rate and formulated the optimization problem 
to maximize the end-to-end data rate. We solved explicitly 
the single relay placement problem with exponential path-loss 
and total power constraint. The asymptotic performance of 
uniform relay placement was also studied for the total power 
constraint. However, the communication model assumed in 
this paper has certain impractical aspects: (i) all the nodes 
should have global knowledge of the codebooks, channel gains 
and transmit powers, (ii) relay nodes need to be full-duplex, 
(iii) global synchronization at block level is required. In our 
future work, we will consider relay placement problems in 
which one or more of these assumptions is relaxed. 



Appendix A 
System Model and Notation 

A. Proof of Theorem [7] 

Let the source and the destination be located at (0, 0) and 
(L,0). There are N number of relays 1,2, ...,A^, where the 
i-th relay is placed at {on^hi). Now let us consider a new 
placement of the relays where for each i G {1,2, ...,N}, the 
i-th relay is placed at (ai,0). Let us define oq := for the 
source and a^v+i '■= L for the destination. This new placement 
results in a reduced (not necessarily strictly) distance between 
any two nodes in the multi-relay channel since now there is no 
variation along y-direction. Let us also assume that the relays 
are enumerated in the increasing order of their x-coordinate 
values. Now let ^ := {i : 1 < i < iV, > L}, B := {i : 
1 < z < iV,ai < 0} and C := {i : 1 < i < TV, a, e [0,L]}. 
Let a minjjg{o,i,2,...,Ar,Af+i},i#j Wi - aj\. Clearly, a<L 
since source and destination are placed L distance apart from 
each other If A is nonempty, place all relays belonging to 
A uniformly (i.e., with equal distance between two successive 
nodes) in the interval \L ~ ^,L\ - Similarly, if B is nonempty, 
place all relay nodes belonging to B uniformly in the interval 
[O, |]. Thus we obtain a new set of relay locations where 
each relay is placed on the line segment joining the source 
and the destination, and the distance between any two relay 
nodes has been reduced. From Equation ([T]), it is clear that R 
increases in the channel gain between each pair of nodes, for 
any given power allocation scheme. Hence the optimum relay 
locations will be such that no relay will be placed outside the 
line segment joining the source and the destination. 



Appendix B 
Single Relay Node Placement: Node Power 
Constraints 

A. Proof of Theorem [2] 

Observe that in Equation (j6| 50,2+51, 2+2^ (1 — a)go.2<?i,2 
is decreasing in a and 050,1 is increasing in a. Also, 50,2 + 
5i,2 + 2^(1 - a)5o,25i,2U=i = 5o,2 +5i,2, which is greater 
than or equal to 050.1 |q=i if 50,1 < 50,2 + 51,2- This motivates 
us to consider two different cases: 50,1 < 50,2 + 51,2 and 
5o,i > 5o,2 + 51, 2- We find that since 50.1 is decreasing in r 
and 51 2 is increasing in r, 50,1 < 50.2 + 51,2 occurs if and 

only it r > rth, where rth — log ' ^ 

Case I: 50,1 < 50,2 - 
maximizer of min{5o.2 + 



" 5i,2- In this case a = 1 is the 
.91,2 + 2^(1 - a)5o,25i,2, "50,1}- 



So the optimum R in this case becomes C 



go.iP 



Further, 



since 50.1 is decreasing in r, in the r > rth region optimum 
relay location is max{rth,0}. 

Case II: .go,i > 50.2 + 5i,2- This holds for r < rth- Now 
as a increases from to 1, 050,1 increases from to 50.1 
and (50,2 + 51,2 + 2y'(l - a)go,29i,2) decreases from (.go,2 + 
51,2 + 2^50,251,2) to (50,2+51,2)- Moreover, both 050,1 and 
(50,2 + 51.2 + 2-\/(l — 0)50.251,2) are continuous functions of 
o. Hence, by Intermediate Value Theorem, there exists Oth G 
[0, 1] such that .go,2+5i,2+2\/(l - 0) 50,251,2 ^ 050, 1 at o = 
Oth. Moreover, min{.go,2 + 51,2 + 2^(1 - 0)50,251,2, 050,1} 
is maximized at o = Oth G [0, 1]. 



From the condition 50,2+51,2 + 2^ (1 - a)5o,25i,2 = "50,1 
we obtain 



1 

90,1 



90,2 



1 - 



91,2 
90,1 



+ A 91,2 1- 



90,2 
90,1 



(11) 



— , ' H 1 ' =. Hence r •) is increasing in u 

2^fe-/c2/«2 2^fc/u-fcV"^ 

over [uth, 1] if and only if its derivative is nonnegative. Setting 
f'{u) to be nonnegative yields the inequality: 



1/u — k/v? 



> -1 



u 
2k 



(14) 



Now let us consider the case —1 



2k 



< 0, i.e u < 2k. 



In this case /(•) is increasing in u. So for r > tq, where 
tq = — ilog(Mo) = ^^log(2fc), /(•) is decreasing in r. Of 
course this condition matters only if vq < rth, since for r > 
rth R always decreases with r. Thus ;/ rg < rth, then the R 
decreases with r for r G [ro,rth]- 



Now for u > 2k, i.e r < ro, from ( 14 1 we will have /'(it) > 
if and only if 

,2 



u 



u 

k 



So /(•) is increasing in u (decreasing in r) if and only 
if w G where it_ and u+ are the two roots of the 

equality version of the above inequality and are given by 7i_ = 

2fc — \/fc,w+ = 2fc + \/fc- Now let us define ?-+ = — - log(it+). 



u- < uth holds if and only if 3\/fc — ^/k + i < 2. It is easy 
to check that this always holds since k G [0, 1] in our case. 
But we are already in w > Uth region. So /(-) is increasing in 
u (decreasing in r) if and only if u G [uthiU+l- 

Again, u+ > itth holds if and only if ^/k + 4 — 3\/fc < 2, 
which holds for any A; G [0, 1]. So > Uth- Hence r+ < rth- 
Thus /(•) decreases with r for r G [r+,rth]- 

Also — 2k + \fk > 2k — uq. Hence r+ < Tq. 

Thus /(-) increases with r for ?■ G [0, r+] and decreases for 
r > r+. If r+ < 0, /(•) decreases with r in [0, L]. 

Summary: 



So for r < rth we have 



max mm < 5o,2+5i,2+2W (1 - oj5o,25i.2, "50,1 

ae[0,l] ■ ' V ' ■ 

Thus our optimization problem becomes 



max 

re [0, rth 




ath5o,i 



(12) 



Let k := e and u := e . Define Uth = e Since 
u strictly decreases with r, our problem becomes 



max f{u) 



where 




" / V ^ V 

Clearly, f{u) is differentiable w.r.t u over [uth,l]- /'("") 



(13) 



1) We found that R is increasing in r for r < 7 + and 
decreasing for r > r+. Let r* be the optimal position of 
the relay node, x* := and x+ := ^ = — ^ log(u+). 

Then x* — max{2;-)_, 0}, where x^ = —j- log I 2e^^ + 



< 1, i-e, A > 
e^2 < 1, i.e.. 



2) rth > if and only if 
log 2. Also r+ > if and only if 2e^^ 
A > log 4. 

3) For A < log 2, we have r+ < and rth < 0. So in this 
case we will place the relay at r* = > rth and hence 
o* = 1 and R* = C{^). 

4) For log2 < A < log4, r+ < and rth > 0. So 
r* ~ and hence we place the node at r* = < rth- 



So o* ^ ath\r=a = 4e-^(l 
5) For A > log 4, r* = r- 



) and hence R* = 



> 0. So substituting the values 
of r+ in Equation (111, we obtain the expression of o* 



(given in Equation ([3])) and thus of R* too (given in 
Equation Q). 

□ 

Remark 1: Let a;th = 77- Then lim;^_5.oo a^th = | and 
liniA-foo x+ = 

Observation 1: ath increases with r in the region r E 
[0,rth]. 

Proof: In the region r < rth, for any fixed r, 

max„g[oa] min |.go,2 + .91,2 + 2^(1 - a)go, 251,2, affo.ij is 

maximized when go,2+gi,2 + 2^ (1 - a)go,2 gi,2 = affoa- Let 
9(0;, r) = 50,2 + gi,2 + 2,/(l - a)go,2gi,2 - a5o,i- Clearly, 
(/(a, r) is strictly decreasing in a and strictly increasing in 
r, and also continuous in both arguments. Also in the region 
r < ''th, '?(«th(?'), = 0. Let us assume that the operating 
point is {r,ath{r)} for some r < rth- Now consider some r' 
such that r' < rth, r' > r. Then q{r' ,ath{r)) > 0. So we 
must increase a from ath{r) in order to bring q{a, r) back to 
0. Hence ath{r') > ath{r)- ■ 

B. Proof of Theorem [i] 

As we observed in the exponential path- 
loss model (see proof of Theorem |2|i, 

min |5o,2 + ffi,2 + 2^(1 - 0)50,231,2, ago,i} is maximized 
for a = 1 if 50,1 < 50.2 + .91.2- This occurs, assuming the 
power law path-loss model, if r^'' < L^'' + (L — r)^'', i.e, 
if r > rth- Note that rth < f . 

In the region r > rth, 5o,i decreases with r, and hence the 
optimum placement point in the interval [rth,i] is rth- 

In the region r < rth, the problem is similar to the one 
of case II in the proof of Theorem [2] Thus we arrive at the 
following problem: 



max f{x) 



where x := 



Xth 



:= ^ and 



1 



(1 



Now f'{x) > if and only if: 

iii-x)-'^-a 



- 1 > 



l)-ri 



(15) 



(16) 



Since x < ^, rj > 1 and the square roots are well-defined, 
the condition above can be written as fi{x) > f2(x) where 
/i(x) = (x-^+i-l)2(l-(i-l)-'') and f2{x) = (l-x)-"- 

— l)^**. Now it is easy to check that fi{x) decreases from 
+00 to as a; increases from to ^. f2{x) increases from 
to {2^ — 1) as a; increases from to ^. So there is a unique 
p € [0, such that f{x) is increasing in a; for x e [0,p] and 
decreasing for x > p. 



Let us now compare p and Xth- Recall that Xth < 5- and 
observe that .Tth verifies: x^^^^ = 1 + (1 — x^''), i.e., {^)^^ + 
(— 1)^'' = 1. Hence, we have: 



/i(a;th) 



{x. 



-j)+i 

th 



/2(a;th) 



—77 

Kh' 1)^ 



(17) 



Clearly, since Xth < 1 and 77 > 1 we will always have 
/i(a;th) < f2{xtb)- Thus we must have p < Xth. This implies 
that R increases with r in [Q,pL], and decreases with r in 
[pL,L]. So R is maximized by placing the relay node at 
r* = pL, where p is obtained by solving /i(a:) ~ f2{x)- 

For the "modified power law path-loss" model, if 2b < L, 
then 50,1 is constant for r e [0, b] but gi 2 increases with r 
for r e [0, b]. So if pL < b, we can improve R by placing the 
relay at b. Hence x* = max{p, in this case. 

□ 

Appendix C 
Multiple Relay Placement : Sum Power 
Constraint 

A. Proof of Theorem |4] 

We want to maximize R given in Equation ([T} subject 
to the total power constraint, assuming fixed relay locations. 



the 



Let us consider Cy-^J2j=iiJ2Ua ^« 
fc-th term in the argument of min{ - } in Equation ([T]l. 
By the monotonicity of C( ), it is sufficient to consider 
X]j=i(X]i=o h'i.k^J Pi.jY ■ Now since the channel gains are 
multiplicative, we have: 

k j-1 



Thus our optimization problem becomes: 




EE 




max mm .9o / 

/ce{l,2,-,7V+l}' ' 'T'O 
j=l ^ 1=0 

N+l 

s.t J2 ^'J ^ 

J=l 
J-1 

and ^P,,j =7jVj e {1,2,--- ,iV + l} 

4=0 



(18) 



Let us fix 71, 72, • • • , 7JV+1 such that their sum is equal to 
Pt- We observe that Pi^N+i for i e {0, 1, • • • , iV} appear in 
the objective function only once: for k — N + 1 through the 

term ^ X^iLo feo"^^ ^ ■ '^i'lee we have fixed 7Ar+i, we need 
to maximize this term over i-^.^v+i, i G {0, 1, • • • , iV}. So we 



have the following optimization problem: Claim 1: Under optimal allocation of 71,72, • • • ,7Ar+i for 



max ^ 



the optimization problem (24i, si = S2 = ■ ■ • = sw+i- 



— r Proof: The optimization problem ( |24[ ) can be rewritten 

i=o as: 



N 

s.t (19) maxC 

=0 



k 



S.t C<&fe^aj7jVfce {1,2,--- ,iV + l}, 



N+l 



By Cauchy-Schwartz inequality, the optimal objective func- 
tion in this optimization proble m is upper bo unded by 

\J (Eto ^»,JV+i)(E.!Io ^) = ^TJV+iEto^ (using the ^ < 

fact that 50,1 = ^o.i ^ {0, 1, • • • , N}). The upper bound is j=i 

achieved if there exists some c > such that = c 7j > Vj G {1, 2, • • • , + 1} (25) 

'»0,i 

Vi e {0, 1, • • • , iV}. So we have The dual of this linear program is given by: 



P^,N+l^ — Vze{0,i,--- ,7V} 
50, i 



s.t 



mill 

fe=l 



Using the fact that EiLo Pi-N+i = Tw+i, we obtain = 

■ Thus, ' jv+i 

ai^6fc/ife + '^i = ^V;e{l,2,... ,7V + 1}, 



1 



k=l 



P^,N+l = (20) > 0, 

E;=o ^ 



90,1 



Here we have used the fact that /lo.o = 1- Now {Pi,N ■ i = 



0, 1, • • • , (/V-l)} appear only through the sum J2fLo^ vT^' {{1j}i<3<n+i,C) which satisfies: 
and it appears twice: for k — N and k = N + 1. We need to 

maximize this sum subject to the constraint Xlilo^ Pt.N = &Va *=C*Vfce|12---A^ + lj 

7Ar. This optimization can be solved in a similar way as ^ /-^'^Hj ~ 1 1 1 ^ h 

before. Thus by repeatedly using this argument and solving 



Mi > 0,i^i > OV; e {1,2,- •• , Ar + 1} (26) 
Now let us consider a primal feasible solution 



optimization problems similar in nature to ( 19 1, we obtain : 



1 



N+l 

E = Pt (27) 

P.,j-= 1 7,-V0<^<J<(jV + l) (21) 

E;=o ^ Thus we have, 5iai7j = (*' i-e-, 7i = 57^- Agam, 62(017^ + 

Substituting for < * < j < (iV + 1) in ([Ts}, we obtain "^72) - C% which implies: 

the following optimization problem : ^ ^ a ^* — (^* 



50, fc 2^ I 7j I Thus we obtain 72 = — j^)- In general, we can write: 



max mm 

/ce{l,2.--- ,N+ 



N+l 



7^ = - (i^-r^) Vfce{i,2,... ,7V + i} 



with := 0. Now from the condition Ea^i^ 7fc = we 



Let us define bk '■= go.k and aj :— Ei=o l^- Observe that oti'^^in- 
bk is decreasing and is increasing with k. Let us define: 



^JV+l 1/1 1 



fe 



Sfc(7i,72,-- - ,7w+i) ■=bk^a3lj (23) 



EiV+i _l_ / 
fc = l afcU;, 6fc_i^ 

- ^Zili^':'^, ^' i g {1. 2, ■ ■ ■ . ^ + 1} (28) 

With this notation, our optimization problem becomes : t* u u u » j *u » • j, ■ ■ ■ • ; *u 

^ ^ It should be noted that smce bk is nonincreasing in k, the 

max min Sfc(7i, 72, • ' ' 1 7Af+i) primal variables above are nonnegative and satisfies feasibility 

ke{i,2,--- ,N+i} conditions. Again, let us consider a dual feasible solution 

N+l 

S.t. ^ Pt (24) 

3 = 1 



({^*, i/*}i<j<jv+i, ^*) which satisfies: 
vl =0V; e {1,2,-- - ,Ar + l} 

JV+l 

k=l 
N+1 

a,^6fe/i^, + i.;=rV/e{l,2,- 

k=l 

Solving these equations, we obtain: 

Q* _ \ 

y-iV + 1 ( A 



,iV + l}(29) 



1,1 ^) 



■spiV+l j_ / 1 1 



j e {1,2,-- - ,Af + l} (30) 



where 



O-N + 2 



:= 0. Since Ofc is increasing in k, all dual 
variables are feasible. It is easy to check that (* = PtO* , 
which means that there is no duality gap for the chosen primal 
and dual variables. Since the primal is a linear program, the 
solution (7i,72,''' i7Ar+iiC*) is primal optimal. Thus we 
have established the claim, since the primal optimal solution 
satisfies it. ■ 

So let us obtain 71, 72,-'' i7a'+i for which si — S2 — 
■ ■ ■ = sjv+i- Putting Sfc = Sfe_i, we obtain : 



k k-l 

bk aj"fj = bk-i ajjj 



(31) 



Thus we obtain 



Ik 



- bk) 1 



Let dh :- 



bk 

. Hence, 



fc-i 



3I3 



k-l 



(33) 



From this recursive equation we have : 



72 = (i2ai7i 



(34) 



73 = ^3(0171 + 0272) = d3ai(l + 02^2)71 (35) 
and in general for fc > 3, 

7fc = rffcain*^;2 (1 + ajdj)-ii (36) 

Using the fact that 71 + 72 + • • • + 7Ar+i = Pt^ we obtain : 



7i 



1 + d2ai + Ylk=5 dkaiU^^^l{l + a^dj) 



(37) 



Thus if si — ,S2 = • • • = sjv+ii there is a unique allocation 
7i,72,''' ilN+i- So this must be the one maximizing R. 
Hence, optimum 71 is obtained by Equation ( [37| . Then, 
substituting the values of {ofc : k — 0, 1, • • • and 
dk : k — 1, 2, • • • , + 1 in Equations (36i and ( [37] i, we 
obtain the values of 71,72, •' ' ,7iv+i as shown in Theorem 

El 

Now under these optimal values of 71, 72,--- ,7Af+i, all 
terms in the argument of min{- • • } in ([T]l are equal. So 
we can consider the first term alone. Thus we obtain the 
expression for R optimized over power allocation among all 
the nodes for fixed relay locations as : R°p^{yi,y2, ■ ■ ■ tUn) = 

C (^5^^^) = C {^^). Substituting the expression for 71 
from Equation ([8]l, we obtain the achievable rate formula (|9|. 

□ 



B. Proof of Theorem [5] 

Here we want to place the relay node at a distance ri from 
the source to minimize < — — \ — so, 1- go, 2 i , Equation 

1 30,1 90,2(1+30,1) I ^ ^ 

|9|. Hence our optimization problem becomes : 



mm s e 

rielOX] 



-pri 



Writing zi — e^"^^ , the problem becomes : 



mm 

Zie[l,eP-f'] 



zi-l 



^pL 



1 



1 



(32) ^^^^ is a convex optimization problem. Equating the derivative 



of the objective function to zero, we obtain 1 



+1 

(^1+1)" 



0. 

Thus the derivative becomes zero at z[ — Vl + eP^ — 1. Hence 
the objective function is decreasing in zi for zi < z[ and 
increasing in zi > zj. So the minimizer is z'l — max{z[, 1}. 
So the optimum distance of the relay node from the source 
is yi=rl~ max{0, 7-']^}, where r[ = ^ log(vTT~e'^— 1). 

Hence ^ = max{i log (^Ve^ + 1 - ,0}. Now > if 

and only if A > log 3. Hence "j; — for A < log 3 and 

$ = ilog(y^ 



1-1 for A > log 3. 



For A < log 3, the relay is placed at the source. Then 

504 = 1 and 50,2 = .91,2 = e"^. Then P04 = 7i = (by 



Theorem 



4i and R* ^ C 



V(e^ + l)CT2 J- 



Also 72 



'1 Pt 
-I 2 



Hence Po,2 = Pi, 2 = 

for A > log 3, the relay is placed at r'l. Substituting the 
value of r'l into Equation (|7ji, we obtain Po.i = 7i = 



Po,2 — 
R* = 

C 



So in this case 



^ c(^SSL^y Since Po,i = ^, we have R* 

1 Pt 



□ 



10 1 be 



C. Proof of Theorem [6| 

For the iV-relay problem, let the minimizer in 
z^,Z2,--- iZ*j^ and let = Mogz^. Clearly, there exists 
i £ {0, 1, • • • , N} such that y*^^ > y*. Let us insert a new 
relay at a distance y from the source such that y* < y < y*_^_l■ 



Now we find that we can easily reach Equation (38 i just by 
simple comparison. For example. 



"i+i 



1 + ZI + --- + Z* 1 + ZI + --- + Z* +ey 



< 



l + zl 



+ 4 



1 + Z* + 



S+1 



1 



First i terms in the summations of L.H.S and R.H.S of 



Equation ( 38 i are identical. Also sum of the remaining terms 



in L.H.S is smaller than that of the R.H.S since there is an 
additional in the denominator of each fraction for the L.H.S. 



Hence we can justify Equation (38 1. Now R.H.S is precisely 
the optimum objective function for the A^-relay placement 
problem ( [TO) i. On the other hand, L.H.S is a particular value 
of the objective in ( 10 1, for {N + l)-relay placement problem. 



This clearly implies that by adding one additional relay we 
can strictly improve from R* of the N relay channel. Hence 

R*{N +1) > R*{N). □ 

D. Proof of Theorem [7| 

When N relay nodes are uniformly placed along a line, 
we will have yk = ]^^- Then our formula for achievable 
,yN) for total power constraint becomes 
7(W) -here /(N) = a + Ef.V T+feS- 
with a — e"+i . 
Since a > 1, we have : 



JV 



rate i?pP*(yi, ?/2, 
: Rn = 



f{N) = a + J2- 



N k 
^ ^ ^ 0*^+1 - 1 



fe=l 

N k 



a + (a - 1)' 



fe=i 
.N 



Thus we obtain: 



liininf /(TV) 

/ pL N , fL 

> lim e"+i H r^(e"+i-l)' 

N-^co \ e«+T 
= 1 



On the other hand, noting that a > 1, we can write: 



N k 



f{N) = a + ia-l)'Y.^kTr—, 

k=l 
N fc+l 

k=l 

= a + ia-lfY^i^ + ^k+r^l 

k=l ^ 
^ f 1 



fe=l 



k=\ 



a + N{a-Xf + (a - l)a~i 



1 - a- 



-AT 



1 - a- 



So 



< 



lim sup /(TV) 

N 

lim (e'fe+iV(e"^ 



l)^ + (l-eT^ 



= 2 



Thus we have limsup^ Rj^ < C (^) and liminf^v Rn > 
C ^ ) • Hence the theorem is proved. □ 
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